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We analyze experimental data of nuclear structure-function ratios F^/F^ and Drell-Yan cross 
section ratios for obtaining optimum parton distribution functions (PDFs) in nuclei. Then, uncer- 
tainties of the nuclear PDFs are estimated by the Hessian method. Valence-quark distributions are 
determined by the F% data at large x; however, the small- a; part is not obvious from the data. On 
the other hand, the antiquark distributions are determined well at x ~ 0.01 from the F2 data and at 
x ~ 0.1 by the Drell-Yan data; however, the large- a; behavior is not clear. Gluon distributions cannot 
be fixed by the present data and they have large uncertainties in the whole x region. Parametriza- 
tion results are shown in comparison with the data. We provide a useful code for calculating nuclear 
PDFs at given x and Q 2 . 

PACS numbers: 13.60.Hb, 12.38.-t, 24.85.+p, 25.30.-c 



I. INTRODUCTION 

Parton distribution functions (PDFs) in the nucleon 
have been obtained by analyzing high-energy nucleon re- 
action data . Such an analysis is crucial for calculating 
precise cross sections for finding new physics phenomena. 
These investigations are valuable for clarifying internal 
hadron structure, and the studies ultimately lead to es- 
tablishment of the nonperturbative aspect of quantum 
chromodynamics (QCD). 

It is well known that nuclear parton distribution func- 
tions (NPDFs) are modified from those of the nucleon Q . 
It was first found by the European Muon Collaboration 
(EMC). Now, major features of the cc-dependent modifi- 
cation became clear experimentally. Although a variety 
of data are not still available in comparison with the nu- 
cleonic case, the PDF parametrization could be done for 

the NPDFs HHI1. The first ^ anal y sis for the NPDFs 
was done in Ref. |4| by using a similar technique to the 
polarized PDF analysis of the Asymmetry Analysis Col- 
laboration (AAC) lg. There are also related studies on 
nuclear shadowing [Tj • The word "shadowing" is used for 
nuclear modification at x < 0.1 throughout this paper. 

These NPDF studies are valuable for describing high- 
energy nuclear scattering phenomena 8]. High-energy 
heavy-ion reactions have been investigated for finding a 
quark-gluon plasma signature. Because such a signature 



should be found in a modification of cross sections, the 
NPDFs should be exactly known. In addition, there is a 
strong demand from the neutrino community to have pre- 
cise neutrino-nucleus, typically the oxygen nucleus, cross 
sections for investigating neutrino oscillation phenomena 
accurately 0, 0, [H| ■ These necessities motivated us to 
investigate the NPDF parametrization. 

In addition, it is interesting to find how well the 
NPDFs are determined. There have been studies of PDF 
uncertainties in the nucleon. It was investigated in the 
unpolarized PDFs and then the studies were ex- 

tended to the polarized PDF uncertainties [dH H3 | . Al- 
though error bands are shown for the NPDFs in Ref. Q , 
they are not based on a rigorous error analysis. Here, 
we calculate the NPDF uncertainties by using the Hes- 
sian method, which is a standard statistical procedure 
for estimating errors 0, 0, Q . 

The purpose of this paper is to report investigations 
after the publication in Ref. 4]. In particular, the fol- 
lowings are added to the previous analysis: (1) Drell- 
Yan data are included in the data set. (2) HERMES 
data are also added. (3) Charm-quark distribution is in- 
cluded. (4) Uncertainties of the NPDFs are estimated by 
the Hessian method. 

This paper consists of the following. In Sec.[H] the x 2 
analysis method, in particular the parametrization form 
and experimental data, is explained. Analysis results are 
shown in Sec. lIIII and they are summarized in Sec. II VI 
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II. ANALYSIS METHOD 

We discuss the \ 2 analysis method. First, the x and A 
dependence of the initial PDFs is explained, and com- 
ments are given on charm-quark distributions. Then, 
experimental data are introduced, and the uncertainty 
estimation method is explained. 



A. Parametrization 



C. Experimental data 



The NPDFs are provided by a number of parameters 
at a fixed Q 2 , which is denoted Qo- The NPDFs could be 
directly expressed by a functional form with parameters, 
which are obtained by a x 2 analysis. However, experi- 
mental data are not sufficient for fixing detailed NPDFs. 
Therefore, it is more practical at this stage to parametrize 
nuclear modification rather than the NPDFs themselves. 
Namely, a NPDF is taken as the corresponding nucleonic 
PDF multiplied by a weight function Wi~. 



f l A (x,Q 2 )=w i (x,A,Z)f i (x,Q 2 ). 



(1) 



The nuclear modification part Wi is obtained by a x 2 
analysis. Here, A is the mass number and Z is the atomic 
number of a nucleus. 

One of the essential points of the x 2 analysis is how to 
choose the x and A dependent functional form. Because 
nuclear modification mechanisms arc different depending 
on the x region, the A dependence could be different in 
each x region. If we would like to describe Wi precisely, 
it could be a complicated function of mixed x and A. 
However, instead of assuming a complicated functional 
form, we use a simple one at this stage. We leave such a 
complicated analysis for our future work. In Ref. jj, a 
simple overall 1/A 1 / 3 dependence is assumed [l^: Wi = 
1 + (1 — 1/A 1 ^ 3 )(x dependent function). Here, we assume 
the same functional form. The weight function used for 
the following analysis is given by: 



w t (x,A,Z) = 1+1- 



1 

A a 



a,i(A, Z) + biX + CiX 2 + diX 3 
(1 - x)P* 



(2) 

where i indicates the parton distribution type, and it is 
taken as i = u v , d v , q, and g. Among these parameters, 
three parameters can be fixed by barvon- number, charge, 
and momentum conservations |J, llq . The motivation is 
explained for choosing this functional form in Ref. £j . 



B. Charm-quark distributions 

In the previous analysis, the flavor number is limited 
to three. However, charm-quark distributions are im- 
portant for practical applications. For example, charmo- 
nium productions are used for searching a quark-gluon 
plasma signature in heavy-ion reactions. The charm dis- 
tributions are also important in neutrino reactions |ll| . 
Therefore, we add nuclear charm-quark distributions into 
the analysis. 

At Q 2 — m 2 , where m c is the charm-quark mass, 
the running coupling constants for the flavor-number 
three and four should agree each other: a s f ~ 3 (m 2 ) — 
a s f ~ (m 2 ). In the leading order (LO), it leads to the 
relation between scale parameters: A3 = A^rric/ A4) 2 / 27 . 
Since the initial distributions in Eq. Q are provided at 
Q 2 which is smaller than m 2 in our analysis, optimized 
parameters for the charm distributions do not exist. The 
distributions appear simply as Q 2 evolution effects. 



In the previous version [4J , the used experimental data 
are limited to the ratios F^/F^ where D indicates the 
deuteron. The data are from European Muon Collabo- 

TABLE I: Nuclear species, experiments, references, and the 
number of data points are listed for the used data with Q 2 > 1 
GeV 2 . 



nucleus 


experiment 


reference # of data 


4 We /n 
ne / yj 


<3T AC PI 
NMP 


I23| 
126] 


16 

1 / 


T i /n 

Be/D 


IN lvlV_^-yO 

SLAC-E139 


f26] 
I23| 


17 

1 T 


P /Pi 


T?MP 88 
PMP QD 

SLAC-E139 


I17| 
[18] 
[23] 


y 
5 
7 




NMC-95 


120] 


1 * 




FN AT F.fifi^ 


rod 
128] 


c 
O 


n /n 

IN / 1 ) 


D l_J 1V1 - O O 

HERMES-03 


[241 
[291 


() 

153 


Al/D 


SLAC-E49 
SLAC-E139 


[21] 
[23] 


18 
17 


Ca/D 


EMC-90 


[18] 


5 




NMC-95 


[26] 


16 




SLAC-E139 


[231 


7 




FNAL-E665-95 


[28] 


5 


Fe/D 


SLAC-E87 


[20] 


14 




SLAC-E140 


[22] 


10 




SLAC-E139 


[23] 


23 




BCDMS-87 


[25] 


10 


Cu/D 
Kr/D 


EMC-93 
HERMES-03 


[19] 
[29] 


19 
144 


Ag/D 


SLAC-E139 


[23] 


7 


Sn/D 
Xe/D 
Au/D 


EMC-88 

FNAL-E665-92 

SLAC-E140 


[17] 
[27] 
[22] 


8 
5 
1 




SLAC-E139 


[23] 


18 


Pb/D 


FNAL-E665-95 


m 


5 


F2/F2 total 






606 


a T7 

(F 2 A /F 2 A ) 
Be/C 
Al/C 


NMC-96 
NMC-96 


m 

[30] 


15 
15 


Ca/C 


NMC-95 


[26] 


24 




NMC-96 


[30] 


15 


Fe/C 


NMC-96 


[30] 


15 


Sn/C 
Pb/C 


NMC-96 
NMC-96 


[3JJ 
[30] 


146 
15 


C/Li 
Ca/Li 


NMC-95 
NMC-95 


[26] 
[26] 


24 
24 


Ff/Ff total 






293 


C/D 


FNAL-E772-90 


[32] 


9 


Ca/D 


FNAL-E772-90 


[32] 


9 


Fe/D 


FNAL-E772-90 


[32] 


9 


W/D 
Fe/Be 
W/Be 


FNAL-E772-90 

FN AL-E866 /NuSea-99 

FNAL-E866/NuSea-99 


[32] 
[33] 
[33] 


9 

8 
8 


Drell-Yan total 






52 


total 






951 



2 



ration (EMC) [H H GjjLthe SLAC-E49, E87, E139, 
and E140 Collaborations^ EJ HI HI, the Bologna- 
CERN-Dubna-Munich-Saclay (BCDMS) Collaboration 
[Hi!, the New Muon Collaboration (NMC) [H, and 
the Fermilab-E665 Collaboration 27, 28]. These data are 
listed in the F^/Ff section of Table[Q 

In addition to these data, we added HERMES data 
for the ratios F A /F^ , where the nucleus A is for nitro- 
gen and krypton [2^. Furthermore, the ratios F A / F A 
{A' ^ D) were measured by the NMC HEESEj, and 
these data are also added. The Drell-Yan data taken by 
the Fermilab-E772 [13 and E866/NuSea |33| collabora- 
tions are added into the data set for the \ 2 analysis. In 
Refs. jsEH^I, Q 2 (dimuon mass) values are not listed. 
Therefore, we calculated the values in the following way 
|34|. Relations between the dimuon mass and the tar- 
get momentum fraction x 2 are listed in Ref. [35|- We 
interpolated these values to obtain the Q 2 information. 

One may note that HERMES 3 He data are not in- 
cluded into the data set. The data are not well repro- 
duced by the present fit, so that the data produce a sig- 
nificantly large \ 2 value. It comes from the fact that the 
3 He is a tightly bound nucleus which cannot be expressed 
by the simple 1 — 1 / A a dependence. In order to reproduce 
such a nucleus, more complicated A dependent function 
should be used for the analysis. 

The number of data points is listed in Table [I] The 
data are for the nuclei: deuteron (D), helium-4 ( 4 He), 
lithium (Li), beryllium (Be), carbon (C), nitrogen (N), 
aluminum (Al), calcium (Ca), iron (Fe), copper (Cu), 
krypton (Kr), silver (Ag), tin (Sn), xenon (Xe), tung- 
sten (W), gold (Au), and lead (Pb). The numbers of the 
F 2l F 2 ( A> ^ D ), and Drell-Yan data are 606, 
293, and 52, respectively. The total number is 951. 

The kinematical range of the used data is shown in 
Fig. □ The smallest x value with Q 2 >1 GeV 2 is 0.0055 
at this stage, and it is rather limited in comparison with 
the proton data (x mm ~ 10~ 4 ) at HERA. The SLAC 
data are taken in the large x, small Q 2 region, and the 
CERN-EMC, NMC, and Fermilab-E665 data are taken in 
the wide x region from small x to large x. The Drell-Yan 
data are in the large Q 2 region. 
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FIG. 1: (Color online) Kinematical range is shown by x and 
Q 2 values of the used data. 



Nuclear modification of the PDFs is expressed by 
the weight functions Wi. We introduce four types by 
assuming the flavor symmetric antiquark distributions 
(u A = d A — s 



q A ) at Q 2 



2. 



u A (x, Qq) = w Uv (x, A, Z) 
d A (x,Q 2 ) =w dv (x,A,Z) 



Zu v (x,Ql) + Nd 

A : 
Zd v (x,Q 2 ) + Nu v (x,Q 2 ) 



q (x,Q ) =Wq(x,A,Z)q(x,Q ), 
g A (x,Q 2 ) =w g (x,A,Z)g(x,Q 2 ). 



(3) 



In the first two equations, the Z terms indicate the pro- 
ton contributions and the N terms indicate the neu- 
tron ones if there were no nuclear modification and 
isospin symmetry could be applied. Although the an- 
tiquark distributions (u, d, s) in the nucleon are different 
[36j, there is no clear data which indicates the differ- 
ence in nuclei at this stage. Therefore, the flavor sym- 
metric antiquark distributions are assumed. The initial 
scale is chosen Qg=l GeV 2 . The MRST01-LO (Martin- 
Roberts-Stirling-Thornc, leading-order version of 2001) 
parametrization j3]} is used for the PDFs, so that scale 
parameter A and charm-quark mass m c are the MRST01 
values in the following analysis. 

Using these NPDFs, we calculate the structure- 
function ratios F A /F A and the Drell-Yan cross sec- 



tion 



ratios <7£>Y/ a DY in the leading order (LO) of 
a s . The NPDFs are given at Ql in Eq. Q, 
so that they are evolved to the experimental Q 2 
points by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution equations in order to calculate these 
ratios. The total \ 2 is defined by 



2l l? 

I° P DY 



* 2 = E 

3 



(4) 



where Rj indicates that the ratios, F A /F A and 
°£>Y'/ jDy- The experimental errors are calculated from 



sys\2 



systematic and statistical errors by ((j^ ata ) 2 = {a a . 
(af at ) 2 . The optimization of the NPDFs is done by the 
CERN subroutine MINUIT [H. 



E. Uncertainty of nuclear PDFs 

Because the situation of the NPDFs is not as good as 
the one of the PDFs in the nucleon, it is especially im- 
portant to show the reliability of obtained NPDFs. The 
uncertainties are shown in the previous version 4]; how- 
ever, they are simply estimated by shifting each param- 
eter by the amount of the error. Of course, a standard 
error analysis is needed for the NPDFs by taking into 
account correlations among the parameter errors. 



3 



TABLE II: Parameters obtained by the analysis. The parameters a Uv , a,d v , and a g are fixed by three conservations. Because 
they depend on nuclear species, they are explained separately in Appendix lAl 



distribution 


a 


b 


c 


d 


u A d A 


fixed (Appendix) 


2.894 ± 0.395 


-9.390 ± 1.068 


7.308 ± 0.866 




-0.3794 ± 0.0461 


8.626 ± 1.551 


-56.64 ± 11.84 


94.11 ± 27.57 


9 A 


fixed (Appendix) 


2.165 ± 3.126 


0.000 (fixed) 


1.349 ± 44.56 



One of the popular ways is to use the Hessian method. 
In fact, it is used for the unpolarized PDF analysis of 
the nucleon [lj] and also for the polarized PDFs [Tj, 13 • 
Because the method is discussed in Ref. 0], we explain 
only a brief outline. 

The parameters of the initial NPDFs in Eq. are 
denoted £j (i=l, 2, • • •, AT), where N is the number of 
the parameters. The \ 2 could be expanded around the 
minimum point £: 

A X 2 = x 2 d + so - x 2 (i) = J2 ii-a-<% > ( 5 ) 

id 



A. Comparison with x-dependent data 

In the actual fit, the parameters for the Fermi- motion 
part are fixed at /3„=/3q=/3 g =0.1 because of the lack of 
large-x data. The parameter a is also fixed at a = 1/3 
Q for the A dependence. The parameters obtained by 
the x 2 analysis are shown in Table ITU Three parameters 
are fixed by the charge, baryon-number, and momentum 
conservations, and they are chosen a Uv , , and a g in 
the analysis. Because these constants depend on nuclear 
species, they are listed separately in Appendix^] An- 
other parameter c g is also fixed since the gluon parame- 



where Hij is called Hessian. The details are discussed 
elsewhere for the uncertainty estimation of the PDFs by 
the Hessian method. For the detailed explanation, one 
may read Refs. 0, 0] about A% 2 and the Hessian. 
A confidence region is identified by providing the Ax 2 
value, which is determined in the following way. The con- 
fidence level P could be chosen as the one-u-error range 
of the normal distribution (P = 0.6826). For one param- 
eter, P — 0.6826 is obtained with A% 2 =1. However, a 
different value should be assigned for the N degrees of 
freedom For example, if there are nine parameters, 
the Ax 2 value is calculated as Ax 2 = 10.427. 

The uncertainty of a NPDF f A {x,^) is calculated by 
the Hessian matrix, which is obtained by running the 
MINUIT subroutine, and derivatives of the distribution: 

M f A M] 2 A JV ( 9f A (x,j) \ tT _ 1 ( df A (x,j) \ 

[5f {x)] =Ax ^[-^) Hij 

(6) 

The derivatives are calculated analytically at the initial 
scale Qq, and then they are evolved to certain Q 2 by the 
DGLAP evolution equations. 



III. RESULTS 

Analysis results are discussed. First, optimized param- 
eters are shown, and x 2 contributions from nuclear data 
sets are listed. Then, fit results are compared with exper- 
imental data. The actual NPDFs and their uncertainties 
are shown for some nuclei at Q\. 



TABLE III: Each x 2 contribution. 


nucleus 


# of data 


x 2 


4 He/D 


35 


56.0 


Li/D 


17 


88.7 


Be/D 


17 


44.1 


C/D 


43 


130.8 
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35 
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57 
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19 


11.8 
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7 
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19 
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24 
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9.8 


Ca/D 
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W/D 
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18.3 


Fe/Be 


8 


6.5 


W/Be 


8 


29.6 


Drell-Yan total 


52 


79.6 


total 


951 


1489.8 
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FIG. 4: (Color online) Comparison with Drell-Yan data of 
/a p D A Y . The ratios (R exp - Rtheoy R th eo are shown 
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FIG. 2: (Color online) Comparison with experimental ratios 
R = F A /F® . The ordinate indicates the fractional differences 
between experimental data and theoretical values: (R exp — 
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FIG. 3: (Color online) Comparison with experimental data of 
i? = F 2 A /F 2 C ' L \ The ratios (fl CTJJ - Rtheoy R theo are shown 

ters cannot be determined easily by the present data. 

The x 2 analysis results are shown in comparison with 
the data. First, x 2 values are listed for each nuclear 
data set in Table ITTT1 The total x 2 divided by the degree 
of freedom is 1.58. Comparison with the actual data is 
shown in Figs. □ GD and g| for the Ff/Ff, F£/F°' L \ 

and Drell-Yan {p V r>Y I ^dy) data, respectively. These ra- 
tios are denoted R ex P for the experimental data and R the ° 
for the parametrization calculations. The deviation ra- 
tios (R exp - R theo )/ R theo are shown in these figures. The 
NPDFs are evolved to the experimental Q 2 points, then 
the ratios (R exp - R theo )/ R theo are calculated. 
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FIG. 5: (Color online) Parametrization results are compared 
with the data of F2 ratios Ff a /F-P and Drell-Yan ratios 



° V dy l aV DY- The theoretical curves and uncertainties are cal- 
culated at Q 2 =5 GeV 2 for the F 2 ratios and at Q 2 =50 GeV 2 
for the Drell-Yan ratios. 

As examples, actual data are compared with the 
parametrization results in Fig. [21 for the ratios F^ a /F^ ) 
and ct^y /& P dy- The shaded areas indicate the ranges of 
NPDF uncertainties, which are calculated at Q 2 =5 GeV 2 
for the F 2 ratios and at Q 2 =50 GeV 2 for the Drell-Yan 
ratios. The experimental data are well reproduced by the 
parametrization, and the the data errors agree roughly 
with the uncertainty bands. We should note that the 
parametrization curves and the uncertainties are calcu- 
lated at at Q 2 =5 and 50 GeV 2 , whereas the data are 
taken at various Q 2 points. In Fig. [S] the smallest- 
x data at ir=0.0062 for F§ a /F® seems to deviate from 
the parametrization curve. However, the deviation comes 
simply from a Q 2 difference. In fact, if the theoretical ra- 
tio is estimated at the experimental Q 2 point, the data 
point agrees with the parametrization as shown in Fig. 

H 



5 



In general, the figures indicate a good fit to the data, 
which suggests that the x 2 analysis should be successful. 
However, there are some deviations as indicated in the 
table and figures. The % 2 contributions are large from 
small nuclei. For example, the Li/D ratios have the x 2 
value 88.7 for only 17 data points. In fact, the Li/D 
ratios in the region x ~ 0.01 deviate from the theoreti- 
cal curve in Fig. The Li/D ratios are measured with 
small errors so that they produce large x 2 values. How- 
ever, if we wish to reproduce the Li/D ratios, the A He/D 
and Be/D ratios cannot be well explained. This is why 
the MINUIT subroutine produced the optimum point al- 
though theoretical calculations deviate from the experi- 
mental Li/D ratios. We also notice that the Sn/C, C/Li, 
and Ca/Li ratios are not well reproduced in the region 
x < 0.04. On the other hand, the figures indicate that 
medium- and large-size nuclei are well explained by the 
parametrization model. 

The Drell-Yan data are taken mainly in the range 
0.02 < x < 0.2 as shown in Fig. H The Drell-Yan 
cross section ratio (J^yl^Y 1S a l mos t identical to the 
antiquark ratio q A (x2)/q A (#2) in the x region, x < 0.1. 
Therefore, the Drell-Yan data are especially valuable for 
determining the antiquark modification in the x region, 
x ~ 0.1. In the smaller x region, the antiquark shad- 
owing is fixed by the F 2 data in any case. Except for 
the W/Be Drell-Yan ratios in the region x ~ 0.02, the 
data are well explained by the parametrization. From 
the constraints of these Drell-Yan cross sections, F 2 shad- 
owing, and momentum conservation, the antiquark dis- 
tributions are relatively well determined in the region 
0.006 < x < 0.1. However, the behavior of the medium- 
and large- a; regions is not obvious. 



B. Comparison with Q 2 -dependent data 

The analysis results are compared with Q 2 depen- 
dent data in Figs. MU\ and for the ratios, F 2 Kr /F 2 D , 
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FIG. 6: (Color online) Q 2 dependence of Ff r ' /F? . The 
curves indicate parametrization results. 
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FIG. 8: (Color online) Q 2 dependence of F 2 S "/F 2 C . 



F 2 N /F 2 , and F 2 S "/F 2 C , respectively. The fit results are 
shown by the curves in these figures. The data are well 
reproduced by the fit except for the Sn/C ratios at small 
and medium x. The tin shadowing is underestimated in 
comparison with the carbon shadowing as indicated in 
the previous subsection. However, we notice that the ex- 
perimental data are not "consistent" in the sense that the 
F 2 Kr /F 2 D and F 2 N /F 2 D ratios tend to decrease at £=0.035 
and 0.045 with increasing Q 2 , whereas the F 2 n /F 2 ratio 
increases. Obviously, more detailed experimental investi- 
gations should be done for clarifying the Q 2 dependence. 
It is especially important for fixing the gluon distribu- 
tions in nuclei. The Q 2 dependence is related partially to 
the nuclear gluon distributions through the Q 2 evolution 



G 



equations. If the experimental Q 2 dependence becomes 
clear, we should be able to pin down the nuclear gluon 
modification. 



C. Optimum parton distribution functions 

We show the nuclear parton distribution functions ob- 
tained by the \ 2 analysis. As a typical medium-sized 
nucleus, the calcium is selected for showing the distribu- 
tions. Because it is an isoscalar nucleus, the and 
are identical. Therefore, u^ a (=d^ a ), q Ca , and g Ca and 
their weight functions are shown in Fig. [5] at Q 2 . 

The valence-quark modification w Uv is precisely deter- 
mined by the data in the medium- and large- a; regions. 
However, the uncertainty band becomes larger in the re- 
gion x < 0.03 although it is constrained somewhat by the 
charge and baryon-number conservations. Obviously, we 
should wait for NuMI (Neutrinos at the Main Injector) 
|39j | and neutrino-factory |4flj projects for clarifying the 
valence-quark shadowing by the structure function F3. 
Although the uncertainties of the nuclear modification 
iy„° are relatively large at x < 0.03, it is not so obvi- 
ous in the valence-quark distribution (xu^ a ), as shown 
on the right-hand side of Fig. because the distribution 
is small in the small- a; region. 

We should mention the possibility that the uncertain- 
ties could be underestimated because we fixed some pa- 
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FIG. 9: (Color online) The weight functions and the nuclear 
parton distribution functions are shown for the calcium nu- 
cleus at Qq. The uncertainties are shown by the bands. 



rameters such as a and f3 in the analysis. In addition, 
there should be uncertainties from the assumed func- 
tional form. These additional factors will be investigated 
in future. In this respect, it is certainly worth while in- 
vestigating the F3 shadowing at future neutrino facilities 
[3^. Eof in spite of the analysis result for the valence- 
quark shadowing in Fig. El 

The uncertainties of the antiquark modification 
are small in the region x < 0.1 because it is fixed by the 
F2 and Drell-Yan data. However, it has large uncertain- 
ties in the x region, x > 0.2. The antiquark distribution 
xq Ca itself is small at x > 0.2, so that it becomes diffi- 
cult to take accurate data for the nuclear modification. 
In order to determine the distribution in this region, we 
need another Drell-Yan experiment which is intended es- 
pecially for large- a; physics |4l| ■ 

The gluon distribution is especially difficult to be de- 
termined by the present data. It is clearly shown in Fig. 
[5] that the modification w^ a and the distribution xg 
have large uncertainties. As explained in the previous 
subsection, the nuclear Q 2 dependence is not clear from 
the data. This fact makes it difficult to fix the nuclear 
gluon distributions. However, we notice that the gluon 
distribution seems to be shadowed although the uncer- 
tainties are large at x < 0.1. 

We notice that the functional form of the gluon weight 
function w g is different from those of the valence-quark 
and antiquark functions, w v and Wq. A similar func- 
tional form was also tested in the analysis. We provided 
a weight function w g , which has the same functional form 
with Wg, as the initial one for the \ 2 analysis without 
fixing the parameter c g . However, the analysis ended up 
with gluon distributions which are similar to the one in 
Fig. [5] It is simply because of the lack of data which are 
sensitive to the gluon distributions. It is the reason why 
we decided to fix the parameter c g in the current anal- 
ysis. The gluon distributions play an important role in 
many aspects of high-energy heavy-ion collisions, so that 
they should be determined by future experimental data. 




FIG. 10: (Color online) The weight functions are shown for 
the nuclei, 4 He, Ca, and Au, at Qq- 
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For example, the eRHIC project [42| could be a promis- 
ing one for determining the nuclear PDFs at small x. In 
order to illustrate the nuclear dependence of the PDFs, 
we show the weight functions for the nuclei, 4 He, Ca, and 
Au, in Fig. Un| 

For general users, a computer code is available on the 
web site i43f f° r calculating the parton distribution func- 
tions for nuclei at given x and Q 2 . The details are ex- 
plained in Appendix B. 



IV. SUMMARY 

The nuclear parton distribution functions and their un- 
certainties are determined by analyzing the experimental 
data of F2 and Drell-Yan data. The uncertainties are 
estimated by the Hessian method. The valence-quark 
distributions are well determined except for the region 
x < 0.03. The antiquark distributions have small un- 
certainties at x < 0.1; however, they cannot be fixed in 
the region, x > 0.2. The gluon distributions have large 
uncertainties in the whole- x region. Obviously, we need 
much accurate scaling violation data or other ones for 
fixing the gluon distributions in nuclei. 
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APPENDIX A: NUCLEAR DEPENDENT 
PARAMETERS 

In Table [H] the constants a Uv , <Xd v , and aq are not 
listed. These constants are fixed by the three conserva- 
tion equations, so that they depend on the mass number 
A and the atomic number Z. For practical usage, we 
express these constants by eight integral values Ii_8 a s 
explained in Ref. Q: 



a„„ (A, Z) = - 
a dv (A,Z) = - 



Zh + {A- Z)h 
Zh + (A - Z)h ' 
Zh + (A - Z)h 
Zh + (A- Z)h ' 



1 



a g (A,Z) = -- 



a Uv {A,Z)\-h 



A^ h 



+ a dv (A,Z){jI 6 



Values of the integrals are listed in Table llVl from the 
current analysis. Using these values together with Eq. 
(|A1(I , one could calculate the constants a Uv , ad v , and aq 
for any nucleus. Then, it is possible to express the nu- 
clear parton distribution functions analytically at Q 2 for 
a given nucleus together with the MRST01 distributions 
[33 in the nucleon. 



TABLE IV: Values of the eight integrals. 



Integral 


Value 


Integral 


Value 


h 


0.2611 


h 


0.3445 


h 


0.1313 


h 


0.1345 


Is 


2.018 


h 


0.2162 


h 


1.016 


h 


0.3969 



APPENDIX B: PRACTICAL CODE FOR 
CALCULATING NUCLEAR PDFS 



One could calculate nuclear PDFs by using the infor- 
mation provided in Appendix A and in Table [HJ How- 
ever, the distributions should be evolved if one wish to 
obtain them at different Q 2 . For those who are not famil- 
iar with such Q 2 evolution, we prepared a practical code 
for calculating the nuclear PDFs at given x and Q 2 . The 
code could be obtained from the web site in Ref. ji^l • 

Instructions for using the code are provided in the 
package. Only restrictions are the kinematical ranges, 



10- 9 < x < 1, and 1 GeV 2 < 



< 10 s GeV . The 



(Al) 



largest nucleus in the analysis is the lead, so that it is suit- 
able to use the code within the range A < 208. However, 
variations of the NPDFs are rather small from A = 208 
to the nuclear matter, one could possibly use the code 
also for large nuclei with A > 208. In the NPDF library, 
we provide the distributions at very small x as small as 
10~ 9 for those who use them in integrating the distribu- 
tions over the wide range of x. However, one should be 
careful that the distributions are not reliable in the region 
x < 0.006, where no experimental data exists. Further- 
more, there is a possibility that higher-twist effects could 
alter the results in the small- a; region. 

The analysis was made in the region, Q 2 >1 GeV 2 , 
where the perturbative QCD is considered to be appli- 
cable. The obtained NPDFs can be used for high-energy 
nuclear reactions with Q 2 >1 GeV 2 . However, there are 
data which are slightly below this region. For example, 
many long-baseline neutrino data are taken in the smaller 
Q 2 region. A useful parametrization was proposed to de- 
scribe the cross section from the deep inelastic region to 
the resonance one |44|. We could possibly make a simi- 
lar analysis in future for describing lepton-nucleus cross 
sections also in the resonance region. 
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